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Abstract 

The Stieltjes constants 7^ appear in the regular part of the Laurent expan- 
sion of the Riemman and Hurwitz zeta functions. We demonstrate that these 
coefficients may be written as certain summations over mathematical constants 
and specialized hypergeometric functions p F p+ \. This family of results general- 
izes a representation of the Euler constant in terms of a summation over values 
of the trigonometric integrals Si or Ci. The series representations are suitable 
for acceleration. As byproducts, we evaluate certain sine-logarithm integrals 
and present the leading asymptotic form of the particular p F p+ \ functions. 
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Introduction and statement of results 

Recently we developed series representations of the Euler constant 7 and values 
of the Riemann zeta function at integer argument, together with other mathemati- 
cal constants, in terms of summations over the trigonometric integrals Si and Ci [8] 
(Propositions 2 and 3, Corollaries 7 and 9). The present work is in a sense a general- 
ization of those results. We present series representations of the Stieltjes (generalized 
Euler) constants that involve sums over certain generalized hypergeometric functions 
p F p+ i. There are underlying connections to Si and Ci and certain logarithmic inte- 
grals of those functions, and our presentation includes some developments of special 
function theory. 

We let C(s) = C( s > 1) be the Riemann zeta function (TTj, ES, Ell El] , T the Gamma 
function, ip = F'/F be the digamma function (e.g., [1]) with 7 = — the Euler 
constant, ip^ be the polygamma functions PQ, and p F q be the generalized hyperge- 
ometric function [2J. Although we concentrate on the Stieltjes constants 7^ = 7^(1) 
corresponding to the Riemann zeta function, we briefly describe how the approach 
carries over to those for the Hurwitz zeta function. (See the discussion section.) The 
Stieltjes constants 7fc(a) [HI EJ El [7j [171 US ESI EH] arise in the regular part of the 
Laurent expansion of the Hurwitz zeta function ((s,a): 

1 00 ( — l)n 

C(*. a) = - + E { -+ln(a)(s - 1)", (1.1) 

where 70(a) = —ip(a). 
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The Stieltjes constants may be expressed through the limit relation [3] 



7n0) 



TV. N—too 



E 

Lfc=0 



ln"(A; + a) \n n+1 (N + a) 
k + a n + 1 



n > 0. 



;i.2) 



Here, a ^ {0, —1, —2, . . .}. The sequence {7fc(a)}fc> has rapid growth in magnitude 
with k for k large and changes in sign due to both k and a. Subsequences of the 
same sign of arbitrarily long length occur. For an asymptotic expression for these 
constants, even valid for moderate values of k, [16] (Section 2) may be consulted. 
The Stieltjes constants appear in applications including asymptotic analysis whether 
in computer science or high energy physics. 
Proposition 1. (a) 



7i = 2£ 



n=l 



1 + -(7 - 2)7 

2 v ' 24 



7T 2 7T 2 n 2 / 5 n n 

— 3 F 4 (1,1,1; 2,2,2, -j-ttV 



+ ln(27m) ( 7 - 1 + - ln(27rn 



+ 5 -7, 



;i.3) 



and (b) the slightly accelerated form 



7i = 2£ 



n=l 



l + -(7-2) 7 , . 

2 w ; ' 24 32tt 4 n 4 



2^2 



3 F 4 ( 1,1,1; 2,2,2,-;--,/ 



+ ln(27m) ( 7 - 1 + - ln(27m 



+ 



73 
144 



7, 



;i.4) 



and (c) 



°° f <-Y 3 7T 2 7T 2 2 

72 = 1 -2(7 + 71) + 2^ {2(1 -7) + 7 2 + y-^ + 7^-3C(3) 



-21n(27m) 



1 - 7 + -| + - ln(2v™)(7 - 1) + - ln 2 (27ra) 



;i.5) 



2 24 2 1 /w 7 6 
Remark. The summand in (1.5) is 0(n~ 3 ) as n — >■ oo. All parts of the Proposition 
may be further accelerated in their rate of convergence. 

Proposition 2. The Stieltjes constant jj may be expressed as a summation over 
n of mathematical constants, terms ln fc (27m), with k = 1, . . . ,j, and the term 

°° / 5 \ 

7T 2 £ ^ 1, 1, • • • , 1; 2, 2, . . . , 2, -; -7T V . (1.6) 

71=1 ^ Z / 

The following section of the paper contains the proof of the Propositions. Section 
3 contains various supporting and reference Lemmas. Some of these Lemmas present 
results of special function theory and may indeed be of occasional interest in them- 
selves. Certain logarithmic integrals of the Si and Ci functions are considered in the 
Appendix. 

Proof of Propositions 

We let P\{x) = Bx(x—[x]) = x—[x]—l/2 be the first periodic Bernoulli polynomial, 
with {x} — x — [x] the fractional part of x. Being periodic, Pi has the Fourier series 
IU (p- 805) 

_gsin^*_ 
Proposition 1. (a) We take a = 1 in the representation 



a~ s a 1_s r°° Pi(x) 



so that 



,, , 1 1 f°° Pi(x) , 
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and 



C(s) + 



1 _ r°° Pi(x) 
(s- l) 2 " ~Jl x s + l 



dx + s 



j; 



Pi(x) 



X 



s+l 



Inx dx. 



Taking s — > 1 in (2.4) and using (1.1) gives 

Pi Or) 



7i 



-i; 



X/ 



\nx dx H 7. 
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(2.4) 



(2.5) 



Here, we have used the well known integral that also results from (2.3) 



x. 



dx 



x/ 



In 



v j ) 2 U + l + 



-7. 



ar 



(2.6) 



From Lemmas 2, 3, and 6 we have 

1 r 00 sin /tz , 1 . tx 1 k 2 ( 5 k 2 

~ J 1 —T- dx = 1 + - (7 — 2)7 — 777 — 777 3P4 ( 1, 1, 1; 2, 2, 2, -; - — 



In k , , 
+ — (27-2 + ln re). 

Taking re = 27m and using (2.1) and (2.5) so that 

1^1/-°° sin(27rnx) , , 1 
7i = r2^-/, ^ lnxrfx + -- 7 , 



(2.7) 



7r n=1 n 71 



(2-* 



gives the first part of the Proposition from (2.7). 

(b) The summand of (1.3) is 0(n -4 ) as n — » 00. By adding and subtracting 
■§2 S^Li 771 = J^^r = 2§8 we transform the summand to be 0(n -6 ) as n — > 00. 



(c) From (2.4) we obtain 
C(s) - 2 



2 y -^Mnx dx- s 1 



s+l 



1 x 



s+l 



In x cte. 



(2.9) 



Then from (1.1) 

poo p ( T ) 

72=/ -A r 1 (2-lax)hix dx, (2.10) 
Ji x 2 

where by (2.5) the term 

poo P ( T ) 

2/ -^Vlnx dx = 1-2(7 + 71). (2.11) 
Ji x 1 

For the term 

r°° PAx) , P 1^1 r°° sin(27rnx) , 9 

/ ^r^\n 2 xdx= Y- — — n J -\n 2 xdx, 2.12 

Ji x 2 n^nJi x 2 ' V ; 

we apply Lemmas 3, 4, and 6. We omit further details. 

Proposition 2. As may be proved from (2.2) by induction, for integers j > 1 we 
have 

cHs, a) = (-i)v- ± Q (j - fc)' ^j^ri + ^r a " s lnJ a 

/•oo p ( T \ 

+(-l) J 7 , , U In^s + a) \j - s \n{x + a)] <fc. (2.13) 
At a = 1 we have 

C (J) (*) = (f^TTI + ^T W ~ lx ^ - slnx ^ dx - ( 2 - 14 ) 
Then by (1.1) we have 

7,- = / — —\n j ~ 1 x(j — \nx)dx 
Ji x l 

= --£- r^^ln^^-lnx)^. (2.15) 
7r ^ n Ji 
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Herein on the right side, j multiplies a contribution from 7^-1. We then appeal to 
Lemmas 2, 3, and 4 for logarithmic-sine integrals, and the result follows. 



Lemmas 

Herein the cosine integral is defined by 

. , f°° cost , , f z cost - 1 , 
Ci(z) = - o?t = 7 + lnz+/ dt. 

Jz t JO t 

As usual, (w) n = T(w + n)/r(w) denotes the Pochhammer symbol. 
Lemma 1. (Hypergeometric form of the cosine integral) 

z 2 ( 3 z 2 \ 

&{z) = 7 + In z - - 2F3 ( 1, 1; 2, 2, -; — J . 

Proof. This easily follows from the expression 

Lemma 2. (a) For a^O, 
Cifaz 



a 

T 



dz — 7 ln(y/x) + - [In (ay) — In (ax)] 
x z 2 

*/ 2 3F4 (l, 1, 1; 2, 2, 2, ?; - x 2 3 F 4 (l, 1, 1; 2, 2, 2, |; 



and (b) for 6 > a > 



/ — In x dx = k {Ci( Kb) (1 + In 6) — Ci(/ta) (1 + In a) — 7 ln(6/ a) 

Ja x 2 



+i[ln 2 M-ln 2 (^)] + ^ 



(3.1) 



(3.2) 



(3.3) 



(3.4) 



' / 3 K 2 b 2 \ o / 3 K z a 

6 2 3F4 I 1,1, 15 2,2,2,-;-— -a 2 3 F 4 ll,l,l;2,2,2,-;- — 



sin/ta,, , . sin/to,, , 
+ 1 + lna) — (1 + ln&). 

a 



(3.5) 



Proof, (a) We first note that l/(£ + 1) = (l)*/(2)/ and (2£ + 2)! = T(2£ + 
3) = 2 • 4 £ (3/2)£(£ + 1)£\. The latter relation follows from the duplication formula 
r(2j + 1) = 4*T(j + l/2)j!/r(l/2). Then by using (3.3) we have 
rv Ci(az), rv dz ~ (-l) £ a 2 ^ 



^ (2£) 2 (2£)! 



(3.6) 



The sums are then rewritten according to 



- ■ ^- 2l 
§ (2£) 2 (2£)! 



(-l)V , 

t7 



1 OO 

4e 



(-i)v 



2£+2 



,2^+2 



4£S(*+l) 2 (2* + 2)!' 



.2 oo 



2£+2 



(2) 2 4' (3/2), (2) £ €! 



= — g-r :;/ • ^1,1,1; 2, 2. 2 - 
(b) We first have, by integrating by parts, 



3 aV^ 



(3.7) 



/ 

Jo 



sin «x 

x 2 



In a; dx 



s 

J a 



sin 

x 



+ k cos /tx In x 



dx sin rex 

X X 



hire 



r b cosK,x, , sin/ta. sin/tfr. 

re/ — — lux cfo + re[Ci(re&) -Ci(rea)]H (1 + lna) — (1 + mfc). (3.8) 



X 



Next, also by integration by parts, 



JO 



cos rex 



lna; cfe = — 



Ci(rea;) 



a X 



dx + Ci(refr) In b — Ci(rea) In a. 



(3.9) 



We then apply part (a) to (3.9) and combine with (3.8 
Lemma 3. (a) For j > 1 an integer, 

00 sin kx , ; , . f°° In 3 ' -1 1 



L 



x^ 



In- 7 x dx = j 



t 



reCi(rei) + 



sin rei 



(3.10) 



(b) With j = 1, 



00 sin 



lnx dx — —k 



Ci(/ct) 



d* + Ci(«) 



+ sin k. 



(3.11) 



Proof. We interchange a double integral on the basis of Tonelli's theorem, 
f^sinra, „• , . / >00 sin/tx [ x Ya?~ l t 



r°° sin /cx . , . f 00 sin kx /" 
/ 5— In- 7 x dx = j / 5— / 

Jl X A Jl X z Jl 



t 



-dt 



-a 



roc [n J - t 1 
Jl t it 

°° ln J ' -1 * r 



. 00 sin kx , 
at / — dx 



-/cCi(/ct) + 



sin /tt 



(3.12) 



with the aid of an integration by parts, (b) follows easily. 

Part (d) of the next result gives a prescription for successively obtaining certain 
needed log-sine integrals, while part (b) complements Lemma 3. 
Lemma 4. (a) 

cos(6x) 



(b) 



s; 

$m(bx) 



\nx dx 



x 



-j; 



Ci(bx) 



x 



dx, 



Inx dx 



1 x^ 



-bd(bx) + 



sm ox 



x 



dx, 



1; 



cos(ox) 



In- 7 x dx 



x 



Ci(ox) 



In 7 1 x dx, 



x 



(d) Let 



Ci(ox) 



with fj (00) = 0, such that 



ln J x dx 



(3.13) 



(3.14) 



(3.15) 



(3.16) 



j f f^{b)db = r w xdx = g . {b)j (3.17) 

J J \ Ob 



with <7j(0) = gj(oo) = and g (b) = sin b — bCi(b). Then for j > 1 



-db, 



and 



where Cj is a constant. 



9i(b) = -jb j 9 -^db + c,b, 



(3.18) 



(3.19) 



(e) Let for j > 1 

m&) = /°° /(^) ^ xdx= ir f{y) w (i) ^ (3 - 2o) 

where /(oo) = and / — >■ sufficiently fast at infinity for the integral to converge. 
Then 

hj(b) = -\ j 'h^(b)db. (3.21) 

Proof. For (a) and (c) we integrate by parts, using that Ci(x) is 0(1/ x) as x — > oo: 
cos(6x) 



ln J x dx 



x 



-'jf 



1 ln i_1 x rfx + Ci(fcx) \n j x 



x 



oo 

1 



X 



(3.22) 



For (b) we use (a) and 



6 Ci(te)^ = 6Ci(te)-^M. 

37 



(3.23) 



so that 



i: 



sin(6x) 



lnx dx 



ar 



s:s: 



cos(fex) 



lnx dxdb 



x 



•b roc 



Ci(bx) 
/o ./> x 

10 



dxdb. 



(3.24) 



For (d) we use the property 



9 W (±) = -iht-i (?) . (3.25) 



db 



Then for j > 1 



6 9 /■« C%) ,__ 7 . 

y 

from which (3.15) follows. For (3.16), we have 



"IzfM ( 3 - 26 ) 



^)= & r^ inj (f)^ (3 - 27) 

from which follows 

J^#) = \ gj {b) - 3 - h9j -M- (3-28) 
The solution of this linear differential equation gives (3.19). That Qj{oo) = follows 
from the Riemann-Lebesgue Lemma applied to (3.17). For (e) we have hj-i(b) = 
—(b/j)df,hj(b) from which (3.21) follows. 

Corollary 1. By iteration of (3.19), we see that gj(b) contains a term with 
l=2£-b\n> +1 b. 

Remark. Imposing #1(00) = gives C\ = 7 2 /2 — 7r 2 /24. 

We note that Lemma 5 just below and the like allows for the straightforward 
integration of the p F p+ i functions arising from the logarithmic-trigonometric integrals 
involved. 
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Lemma 5. For p ^ 0, q ^ — 1, and (q + l)/p ^ — 1, 



J n q 3 F 4 (ai, a 2 , 03; h, b 2 , b 3 , c; -n p /A)dn 

= — — 7 4^5 ai, a 2 , a 3 , ; 61, b 2 , 63, c, h 1; — - . (3.29) 

g + 1 V P P 4 / 

Proof. We have 

/ «• 3F4K a,, a 3 ; ft,, ft 2 , ft,, c; = £ fr^jMffi 1 (-1)' I / K "+V K 



We then use 



V v it 

to complete the Lemma. 

We then call out a special case at c = (q + I) /p. 
Corollary 2. For p ^ 0, 



(3.31) 



/ 3 F 4 (a, a, a; b, b, b, - ; | dn 

J \ p 4 J 



K q+l ( q + 1 k p \ 

= — — 3-^4 a, a, a; b, b, b, hi; — - . (3.32) 

q + 1 V P 4 / 

In particular, this result applies with p = q = 2 when integrating the function of (3.6) 
with respect to a. 

Lemma 6. (Leading asymptotic form of special p F p+1 functions) For z — > 00 we 
have (a) 

2 F 3 (l, 1; 2, 2, 1 -z) ~ - z (7 + In 2 + ± In , (3.33) 
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(b) 



Fa 1,1,1; 2, 2, 2, 



6(1 -7) + 37 2 - y -61n2 + 67ln2 + 31n 2 2 



31n 2 z 



+ -(7- 1 + ln2)lnz + - 

Z 4: Z 



(3.34) 



and (c) 



F R 1,1,1,1:2,2,2,2, 



12(7 - 1) - 67 2 + 2 7 3 + — - |tt 2 -6 In 2 2 + 67 In 2 



3 3 

+21n 3 2 + 61n2(7- 1) \nz + 3 In 2 2 \nz - - In 2 z + -(7 + In 2) In 2 2 + In 3 z 



7T 



+ (6(l-7) + 3 7 2 --)ln(42:)+4C(3) 



(3.35) 



Proof, (a) This follows immediately from the definition of Ci in (3.1) and the 
relation (3.2). However, we use a more general procedure based upon the Barnes 
integral representation of p F q ([20J, Section 2.3). We have 



2-^3 



Fx f l,l;2,2,-;-z 



— / T(-s)T(s + l)g(s)z s ds, (3.36) 

Z7TZ J— zoo 



where the path of integration is a Barnes contour, indented to the left of the origin 
but staying to the right of —1, and 



r(s + l)g(s) = V (-) 



3\r 2 ( s + i) 1 



(3.37) 



T 2 (s + 2) T(s + 3/2) \2j (s + l) 2 T(s + 3/2) ' 
The contour can be thought of as closed in the right half plane, over a semicircle 
of infinite radius. We then move the contour to the left of s = — 1, picking up the 
residue there. We have 

d T(-s) T(-s) r 



dsT(s + 3/2) T(s + 3/2) 
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-s)+V(s + 3/2)], 



(3.38) 



and '0(1/2) = —7 — 2 In 2. Therefore, we have for s near —1 

r (!) F(JT§2) = \ + h + ln2 >< s + : > + °l< s + " 2 1- < 339) 

with z s = l/z + In z(s + l)/z + 0[(s + l) 2 ] and find 

rf-Ves s — 1 7 \{~ S } ZS 1 \ = - (-i + \n2 + -\nz\ , (3.40) 

\2J s ~ (s + l) 2 r(s + 3/2) A' 2 7' V ; 

giving (3.25). 

(b) We proceed similarly, with 

(5 \ i fioc 

1,1,1; 2, 2, 2, -; -zj = — T(-s)T(s + l)g(s)z s ds, (3.41) 

and 

r <* + - r (I) s^u+v 

Again using an expansion like (3.39) and evaluating the residue at s 
(3.34). 

(c) Similarly, we use 

1, 1, 1, 1; 2, 2, 2,2,-;-zj = —J ^ T(-s)T(s + l)<?(s)^s, 

with 

r« + DM - r (|) ± j _ I _ = r (§) (< - ^ _ _ (3 . 4 4) 

and compute 

V27 s " (s+ l) 4 r(s + 3/2) 
14 



(3.42) 
= — 1 gives 

(3.43) 



Remarks. We note that results like part (b) of this Lemma allows us to find the 
limit as y — > oo for the integral of Lemma 2. 

We have presented the algebraic part of the asymptotic form of the p F p+ i func- 
tions, that is the leading portion. The higher order terms come from the exponential 
expansion of these functions, and they are infinite in number. According to ([20J, 
Section 2.3, p. 57) this expansion has the form 

OO A 

E(-z) = (2^ 2 )V^^^^, (3.45) 

where 9 = 1 — p — a and A = (2n)~ 1 / 2 2~ 1 / 2 ~ 6 '. Here, a is either 3/2 (for the p = 2 
case) or otherwise 5/2. The other coefficients come from an inverse factorial 
expansion of the function g(s), 



g(s) = 2-A s 



f> Aj 0(1) 



(3.46) 



, =0 T(2s + 1-9+ j) T(2s + 1-9 + M) 
The full asymptotic expansion of the p F p+ i functions includes the algebraic portion 
H(—z) given in the Lemma, together with (3.42), as 

P F P+ i (l, 1, • • • , 1; 2, 2, . . . , 2, | -z) ~ E{-z) + E(z) + H(-z). (3.47) 

Lemma 7. (Alternative integral representation) 

\ f x 2 \-ncos(e~ x/2 K) + e x/2 sm(e~ x/2 K)]dx. 

K JO 

(3.48) 

Proof. This follows by using the special case Laplace transform 



3 ^(Vl,l;2,2,2,-;-^ 3 ' /X 



5 5 ? i i n : 



T 4 / 2 



xV^dx = (3.49) 
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Then, with the integral being absolutely convergent, we may interchange summation 
and integration in 

3 F 4 2, 2, 2,5; 4) = Igjf- A -^_l_(_£y 1*. (3.50) 

and the result follows. 

We mention the connection between integrals expressible as certain differences 
of the confluent hypergeometric function \F\, differences of the incomplete gamma 
function j(x,y), and a i_F 2 function. We have 
Lemma 8. We have for Re \x > — 1 

/ x 11 ' 1 sm(ax)dx = [ 1F1U1; fj, + 1; ia) — iFi(u; fj, + 1; — ia)\ 

Jo 2/z 

= - jt* )""!'!'^.-' ) - (-1)"7(/M<>)] 



/x+ 1 V 2 '2 2 ' 4 
We note that by logarithmic differentiation with respect to fi this result generates 

a family of logarithmic sine integrals. I.e., for Re \i > — 1, 

/ sin(ax) ln fc x = [ — — ] / x^" 1 sm(ax)dx. (3.52) 

JO \OjJi) Jo 

Proof. The first line of (3.51) is from ([12], p. 420, 3.761.1). The second line 
follows from the relation ([12], p. 1063) 

7(a, x) = — x a iFi(a; a + 1; — x). (3.53) 

a 

For the third line we write 

l F 1 ( j i; 11 + 1; ia) - ^Gu; /i + 1; -ia) = £ (/ ^\ ^ [1 - (-1)*] 
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oo / \ n 2m+l 

= 2z V - ^ hrr ; +l -g -. (3.54) 

^ (/i + l) 2m+1 (2m + l)! 1 ; 

The Lemma is now completed with the use of the duplication formula, as at the 
beginning of the proof of Lemma 2, and the relation 



(»2m+l \l 1 \l \2* 



(fi + l) 2m+1 2(///2 + m + l/2) 2(/i + l)fe N ) 
Discussion 



(3.55) 



We mention how hypergeometric summatory representations may be obtained for 
the constants 7fc(a). When k = we must recover from (2.2) the known representa- 
tions (e.g., [11], p. 107) 

lnT(a) = (a-l-)\na-a + l\n(27i) - f° ^-dt, (4.1) 
V 2 / 2 7ot + a 



V»(o) = -7o(o) = In a - — + / ^ , J^ *' ( 4 - 2 ) 



and 

j_ + Px(t) 
2a Jo (t + a) 

leading to, for Re a > 0, 

oo 

■0(a) = —70(a) = In a — - — h [2 cos(27rja)Ci(27rja) — sin(27rja) [it — 2Si(27rja)]] . 

a j=i 

(4.3) 

Matters are more interesting, and challenging, for k > 1. 
From (2.2) we have 

a~ s lna a 1_s lna a 1_s 



2 8-1 (s-1 



|2 



+ f(^^l sln(a; + a) - 1 ^ (44) 
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so that 



In a In 2 a r°° PAx) ri . . „. , . , . 



(x + a) 

Again the Fourier series (2.1), applies so that the term 

f°° PAx ri . , , 1^1 Z" 00 sin(27rnx) M . . , 

/ . u . J in x + a - 1 da; = V - / — 1 - In (x + a )}dx 

Jo (x + a) 2 7r ^ n jo (x + a) 2 

n^nJa x 2 v y 



1 00 1 r 
n^nJa 



[sin(27rnx) cos 27ma — cos(27rnx) sin 27ina] 



l-lnx)dx. (4.6) 



x* 

A special case occurs at a = 1/2, when ((s, 1/2) = (2 s — l)C(s), implying that 
7i (1/2) = 7 i - 2 7 ln2 - In 2 2. We have from (4.4) and (4.5) 

, /^-.\ , „ 1, o„ 1 ^ (-l) n f°° sm(2imx) , 
- 7l (l/2) = ln2 + -ln 2 2-- ^/ — — ^(1-lnxdx. (4.7) 



n=l 

Here 



(°° Sin{27 l nx) dx = -«Ci(«/2) + 2 sin(«/2), (4.8) 

il/2 X 2 



r°° sin(27rnx) 

'1/2 

with the last term vanishing when k = 2nn. The logarithmic integral in (4.6) can be 
decomposed as J-jy 2 = fl/ 2 + Ji°° and Lemma 2 and (2.7) can be applied respectively. 
This leads to a hypergeometric summatory expression for 71 (1/2). 

Lastly, we mention that our methods for treating logarithmic integrals may be 
applied to other functions, including the Bessel functions of the first kind J n . These 
integrals may be of the form 

^l\n k (bx)dx, 

x J 
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and in particular for n — 1 we expect analogies with our previous trigonometric 
logarithmic integrals. However, we do not pursue that topic here. 

Summary 

We have presented summations over particular p F p+ i generalized hypergeometric 
functions and mathematical constants giving the Stieltjes constants jk- We have de- 
scribed the asymptotic form of the p F p+ i functions with repeated parameters, and 
indeed the extended asymptotics may be used to accelerate the convergence of the 
hypergeometric-based summations. We have considered certain logarithmic-sine inte- 
grals. We showed that these are essentially certain iterated integrals with the cosine 
integral Ci at their base. Our results generalize significantly simpler summatory ex- 
pressions for the Euler constant 7 in terms of evaluations of Ci or the sine integral 
Si. 
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Appendix: Logarithmic integrals 

The following two integrals are given in the Table [12] (p. 647): 



oo n 

Ci[x) Inx dx — —, 



(Al) 



and 



POO 

/ si(x) In x dx — j + 1. 
Jo 



(A2) 



Here, 



si(x) 



sint 



t 

We show that these two integrals follow as special cases of the following. 
Proposition Al. For a > 0, 



(A3) 



CO I 

e~ at cos t(l - In t)dt = — — <| 2 cot -1 a + a 

ii 2(1 + a 2 ) 



ln(l + — ) +2(l + 7 + lna) 



(A4) 



and 



POO 

/ e~ at sin til - In t)dt = — - 
7o V ; 2 1 



-2a cot x a + ln (l + — ) +2(1 + 7 + lna) 



(A5) 



We first observe the special cases under formal interchange of integrations 



POO 

/ Ci(x) 
Jo 



■oo poo cos f 

mx dx = — I mx dtdx 

Jx t 



00 cos t r f 



o t Jo 



lnx dxdt 



POO 

/ cos£(l — hit)dt, 
Jo 



(A6) 
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r°° , s , , f°° , f 00 sint , , 
lnx ax = — / mx dtdx 



lo Jx t 

'- 00 sint /•* 

mx axat 



and similarly 

/•oo 

/ si(x) 

JO 

« = „ _ r°° srnt 

Jo £ JO 

= / sint(l - ]nt)dt. (A.7) 
Jo 

Indeed the integrals (A. 6) and (A.7) are not convergent. However, the a — > limit in 
the Proposition gives the integrals (A.l) and (A. 2). 

Proof. We just prove (A. 4), as (A. 5) goes similarly. From the Gamma function 
integral 

/ e- at t z - l dt = a- z T(z) (A.8) 
Jo 

we have 

/ e- at t z - 1 \ntdt = —a- z T(z) 
Jo oz 

= a~ z T(z)[iP(z) -In a]. (A9) 

We recall the relation ip{n + 1) = H n — 7, where H n is the nth harmonic number, and 
then for integers j > 

roc 

/ e- at t 2] \nt dt = a- 2j+1 T{2] + l)W2j + 1) - In a] 
Jo 

= a- {2j+1) (2j)\[H 2j - 7 - In a}. (A 10) 

Then integrating term-by-term we have 

roo 00 ( 1 \3 roo 

/ e~ at cos tlntdt = Y — — ^— / e- at t 2] In t dt 
Jo (2j)\ Jo 
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= Eyyr«' (2i+1) (2i)![%-7-l.io] 

= _(i±M + g ( _ 1)3//2ja - w+1) . (A11) 

1 + a j=0 

The latter sum may be performed with the substitution of an integral representation 
for the digamma function, or by manipulating the generating function 

gc-i)-^- = -^a±£), ( Ai2) 

n=0 1 "I" Z 

and recalling that cot -1 z = (l/(2i)) ln[(l + iz)/(l — iz)]. We have 

+ = 2 1 = _M1±£> + Mi_£>. (A13) 

n=0 m=0 1 + Z Z I 

We then put z = ia to obtain the sum of (A. 11). Using the elementary integral 
J °° e" at costdt — a/(l + a 2 ) completes the Proposition. 
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